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2 ■ Abstract 
ON 



We study coherent dynamics of tight-binding systems interacting with static 
and oscillating external fields. We consider Bloch oscillations and Wannier- 



^3 

£5 ' Stark localization caused by dc fields, and compare these effects to dynamic 

o " 

localization that occurs in the presence of additional ac fields. Our analysis 
relies on quasienergy eigenstates, which take over the role of the usual Bloch 
waves. The widths of the quasienergy bands depend non-monotonically on 
the field parameters. If there is lattice disorder, the degree of the result- 
ing Anderson localization is determined by the ratio of disorder strength and 
quasienergy band width. Therefore, the localization lengths can be controlled, 
within wide ranges, by adjusting the ac amplitude. Experimental realizations 
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of our model systems are given by semiconductor superlattices in far-infrared 
laser fields, or by ultracold atoms in modulated standing light waves. In both 
cases the system parameters, as well as the field amplitudes and frequen- 
cies, are readily accessible to experimental control, suggesting these as highly 

attractive candidates for systematic study of localization phenomena. 
72.15.Rn, 73.20.Dx, 42.50.Hz 
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both fields is given by g^^ 



I. INTRODUCTION: THE "RENORMALIZATION" OF THE LANDE G-FACTOR 

More than 25 years ago Haroche, Cohen- Tannoudji, Audoin, and Schermann demon- 
strated experimentally that the Lande factor gp of atomic hyperfine levels F can be mod- 
ified by an oscillating magnetic field Working with Hydrogen and Rubidium atoms, 
respectively, and applying an oscillating field B\ cos(a;t) perpendicular to a static field B, 
the authors showed that the effective value g^^ of the Lande factor in the presence of 

gpJo(gpfiBBi/u), where /xg is the Bohr magneton, and J 
is the ordinary Bessel function of order zero. Similar results have also been published by 
Chapman [0] and Yabuzaki et al. [||,|J. 

This renormalization of the Lande factor by oscillating fields is closely related to the 
subject of the present paper: the renormalization of the band structure in tight-binding 
lattices caused by external ac electric fields, and the resulting possibility of manipulating 
tunneling and localization effects by suitably tuning the ac amplitudes. In order to elucidate 
how an oscillating field can affect the Lande factor, and to prepare for the following discussion 
of tight-binding systems, we consider an electron in a state with total angular momentum 
F = 1. We assume that a static magnetic field of strength B is applied in the ^-direction, 
together with an additional static field Bq and an oscillating field B\ cos (cut) in the re- 
direction. The Hamiltonian, restricted to the F — 1 subspace, can then be written as 



'lOO^ 



H(t) = g lf j, B B 




0-1 



gifjLs [B + B 1 cos(ut)] -j= 



^0 10^ 
1 1 



(i) 



Because H{t) is periodic in time with period T = 2n/u, there are Floquet states [§-0] 



ip a (t) = u a (t) exp(-ie a t) 



(2) 



where the T-periodic Floquet functions u a (t) solve the eigenvalue equation 



(H(t) - id t ) u a (t) = e a u a (t) 



(3) 
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The eigenvalues e a are called quasienergies. 

Now there is a simple, but important detail to note. Suppose that Uj(t) is a T-periodic 
solution to (|3|), with quasienergy Ej. Then, for every positive or negative integer m, the 
function Uj(t) exp(imut) is also a T-periodic eigensolution, with quasienergy Ej+muj. Thus, 
the quasienergy spectrum consists of "Brillouin zones" of width u, and the index a in ([3]) 
should be regarded as a double-index: a = (j,m). The eigenvalue problem (|3]) is defined in 
an extended Hilbert space of T-periodic functions ||, with scalar product 

«•!•» = y X^H-)- ( 4 ) 

This is the usual scalar product combined with time-averaging. All eigenf unctions 
Uj(t) exp(imLdt) are needed for the completeness relation in the extended Hilbert space. 

On the other hand, all solutions that differ only by a factor exp(imu)t) fall into the same 
class, in that they represent the same physical state ipj{t) = Uj(t) exp(-iSjt). An arbitrary 
solution ip(t) to a time- dependent Schrodinger equation with a T-periodic Hamiltonian H(t) 
can be expanded according to 

ip(t) = y~] ajUj(t) exp(— ie-jt) , (5) 

3 

where j labels eigensolutions of (|^) that correspond to different states. This expansion shows 
that for quantum systems governed by a time-periodic Hamiltonian the Floquet states and 
quasienergies play the same role as energy eigenstates and energies do in the static case. The 
coefficients cij do not depend on time. Note that the Fourier index m does not appear: only 
one representative Uj(t) of each class of eigenfunctions is required here. For this reason, we 
will not use the double- index notation in the following, and merely indicate the Brillouin- zone 
structure of the quasienergy spectrum by adding "mod u" to the corresponding formulae. 

We now compute the quasienergies for the Hamiltonian (|l|) perturbatively. A unitary 
transformation yields 



4 



H(t) = 9l fx B B 



V2 



^0 10^ 
1 1 



gxlis [B + B x cos(u;£)] 




0-1 



(6) 



We will denote the first term on the right hand side as H^ B \ If B = 0, the wave functions 



o 



V 



exp -igifi B (B Q t + 



CO 



V / 



V 1 / 



( ■ /n Bi$m.Uot). 
exp I +ig 1 /j, B {B t H 
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(7) 



solve the time-dependent Schrodinger equation. Obviously, these are Floquet states with 
quasienergies Em — fl'i/^B-Bo • M for M = —1, 0, +1, to be taken modulo uj. 

If B is small, B <^ B w Bi, the matrix can be considered as a perturbation. 

Perturbation theory for quasienergy eigenstates in the extended Hilbert space works exactly 
like the well-known Rayleigh- Schrodinger perturbation theory for energy eigenstates, pro- 
vided the proper scalar product @ is employed. In our case, the "unperturbed" Floquet 
functions Uj(t) can be identified as the T-periodic parts of the wave functions (^). If there 
is no integer n such that gi/i£-E?o = nu, then there is no resonance, and we can resort to 
nondegenerate perturbation theory. For small B it suffices to calculate the diagonal elements 
of the perturbation. Since 



(( Uj (t)\HW\ Uj (t))) = 



(j = 1,2,3) 



(8) 



there is no first-order correction to the quasienergies. 

The situation is different if giUsBo = nuj for some integer n — i.e., if the energy of n 
photons^ coincides with the energy difference between adjacent Zeeman states. Then the 



1 Since the ac field is treated as a classical external force, rather than a quantized field, the term 
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Brillouin-zone structure of the quasienergy spectrum becomes crucial: the Floquet eigen- 
functions 



«i(t) 



Bi sm(cut)\ 
exp | -igmB , u 2 {t) 



V 



Anuit 



U 3 (t) 



V 



V 1 / 



cxp 



/ B 1 sm(ut) 
\+igifx B 



(9) 



are degenerate. Therefore one has to diagonalize the matrix of the perturbation within 
the degenerate subspace. Using the generating function for the ordinary Bessel functions J k 
of integer order, 



k=~oo 



this matrix becomes (with j, k = 1, 2, 3) 
(( Uj (t)\H^\u k (t))) 



V2 



gifiBB J n 



(giliBBx \ 



Hence, the quasienergies are given by 



■D n J, 



( gifj>BBi 

\ CO 



) gi^sB ■ M mod 



for M = 0, ±1. Thus, in effect the g-factor is "renormalized" by the ac field: 



9\ 



(10) 



(11) 



(12) 



(13) 



Originally, the renormalization of the atomic g-factors had been explained within the 
dressed-atom picture [PJlOfl, which amounts to treating both the atom and the oscillating 



"photon" is, strictly speaking, not justified here. We use it nonetheless to express the resonance 
condition in a form that can be easily remembered. 
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field quantum mechanically. Within the Floquet picture the field is not quantized, but is 
regarded as an external force acting on the atom. The two approaches give identical results, 
as long as the ac fields are sufficiently intense || . The line of reasoning applied to the above 
example is essentially a formalized version of a semiclassical argument given earlier by Pegg 
and Series Jll] . 



The lesson to learn from this example is that resonant oscillating fields can substantially 
alter the level structure of atoms or molecules. Similarly, ac fields can modify the band 
structure in periodic potentials. This effect is of possible relevance for present experiments 



on semiconductor superlattices exposed to far- infrared radiation |1^-|15|: ideally, one can 
alter the quasienergy (mini-)bandwidths in these systems by tuning the amplitude of the 
radiation field, and thereby control tunneling and localization phenomena in the laboratory. 
In particular, experiments of this kind might shed new light on Anderson localization. 

In the following, we explain this possibility in more detail. We first discuss a simple 
two-level problem, which already exhibits many of the central features, and then move on to 
one-dimensional tight-binding lattices. We review the phenomena of Bloch oscillations and 
Wannier-Stark localization occuring in the presence of dc electric fields, and the dynamic 
localization that one finds in combined ac and dc fields. A crucial step then follows: we 
introduce disorder into the tight-binding model, and demonstrate that the parameter which 
determines the degree of disorder-induced localization is the ratio of disorder strength and 
the quasienergy bandwidth. The case of an isolated defect is treated analytically; numerical 
results are given for randomly disordered lattices. The role of interband effects is discussed 
within the scope of a two-band model. The paper then closes with a brief discussion and 
outlook. 

II. THE TWO-LEVEL SYSTEM 

Let us first consider a rudimentary "lattice" with just two sites, each of which supports 
only a single quantum state — i.e., a two- level system, which interacts with an external dc 
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electric field of strength Fq and an additional ac field of strength F\ and frequency uo: 

Htis(t) = jv z + {J,[F + Ft cos(ut))a x . (14) 

The energy difference between the unperturbed states is denoted by A, /i is the dipole matrix 
element, and a x , o z are the usual Pauli matrices. This Hamiltonian describes, for example, 
an electron in a semiconductor double-well heterostructure under dc bias in the presence of 



far- infrared radiation |16| , |17| , provided the dynamics remain restricted to the lowest doublet 
of energy eigenstates. 

Evidently, H t i s (t) is of the same type as the Hamiltonian (|l|) (but with two levels, or 
effective spin 1/2 rather than spin 1), and the perturbative calculation of the quasienergies 
proceeds exactly along the lines of the previous example. For A = the Floquet functions 
take the form 



1 ( 1 ^ A^i 

-= exp[=Fi— siii(wt)] , (15) 



and the quasienergies are 

e± = ±/xF mod u . (16) 

If the ac field is not resonant, i.e., if 2/zF ^ nu for all integer n, these quasienergies remain 
unchanged to first order in A/u. However, if 2fiFo = nuj for some integer n, then degenerate 
perturbation theory in the extended Hilbert space yields 

If there is no dc field, F Q = 0, the resonance condition is fulfilled automatically with n = 0. 
In this case the two quasienergies cross when 2/zFi/u; is approximately equal to a zero of 
the Bessel function J . This fact, which had already been derived by Shirley in 1965 
and thoroughly studied experimentally 0, has sparked renewed interest recently, after it 
was pointed out by Grossmann et al. [|T^-|2D| that tunneling in one-dimensional symmet- 
ric double-well potentials can be coherently suppressed by an oscillating force. Assuming 
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that the wave packet dynamics can be captured by a two-level approximation, an obvious 
necessary condition for the suppression of tunneling is the crossing of the two quasiener- 
gies emerging from the lowest doublet of energy eigenstates of the unperturbed double well, 
and Shirley's formula (|T7|) gives a reasonable approximation to the parameters where the 
quasienergy crossing occurs [[H],^2|]. The two- level approximation becomes insufficient when 
the ac frequency or amplitude is so high that the field not only mixes the lowest two dou- 
ble well eigenstates but also couples these states to higher doublets. Then the quasienergy 
"ground state" splitting can become even larger than the bare tunnel splitting A. This fact 
enables fast, efficient, and externally controllable population transfer in ac-driven double 
wells [p3|l . It should be pointed out, however, that the full quasienergy spectral problem 
for periodically driven anharmonic oscillators is technically far more demanding than the 
computation of spectra for iV-level systems. At present, it is even unknown whether a 
periodically forced oscillator with a quartic anharmonicity has a singular or an absolutely 
continuous quasienergy spectrum, if the amplitude of the driving force is large |24| . 



III. THE INFINITE TIGHT-BINDING LATTICE IN DC FIELDS: 
WANNIER-STARK LOCALIZATION AND BLOCH OSCILLATIONS 

What happens to the approximate spectrum flT7| ) when the number of lattice sites is 
increased from two to infinity? This question leads us to consider the tight-binding Hamil- 
tonian 

ffo = -yE(|* + l><<| + \£>(l+l\) (18) 
with scalar potential interaction 

H mt (t) = e[F + F 1 cos(ut)]d £ \i) t (£\ . (19) 

e 

The system defined by the Hamiltonian Hq + Hi nt (t) can be taken as a model for an electron 
in an ideal superlattice in the presence of both terahertz radiation, linearly polarized along 
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the growth direction of the lattice, and additional dc bias. An "atomic" (Wannier) state at 
the E-ih site is denoted by \£), and d is the lattice constant. If there is no external electric 
field, F = Fx = 0, then the Bloch waves 

Xk = Y, e ~ tkM \Z) (2°) 

i 

solve the stationary Schrodinger equation, H Xk = E(k)xk, with energy eigenvalues 

E(k) = -^-cos(kd) . (21) 

Thus, A is the width of the resulting energy band. 

If there is a dc field F Q , but no ac field, the Wannier-Stark states 

^M^)^ (22) 



constitute a complete set of energy eigenstates Their eigenvalues, 



E m = m ■ eF d , (23) 

form the so-called Wannier-Stark ladder. A particular Wannier-Stark state <p m is localized 
around the m-th site, with a localization length of the order of A/(eFod) lattice periods. 
The addition theorem for the Bessel functions Jk(z) of integer order, 

+oo 

E Mz)J k+P (z)e ika = J p (2zsm(a/2))e ip ^ 2 , (24) 

k=— oo 

directly reflects, for a = 0, completeness and orthogonality of these states. 
The propagator for the dc-biased lattice reads [p5[ 



K i,i'{t) = J2(£\(p m )((p m \£') exp(-iE m t) 

m 

= J '-" (^5 S ' n irr) ) CXP ( !(/ " ef-^f^ - U'eF dt\ , (25) 



where fl24] ) has been used. Because the Wannier-Stark spectrum (|23| ) is a ladder with 
equidistant level spacing AE = eFod, the propagator K^if) is periodic in time with the 
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Bloch period TbiocH = 2tt/ (eF^d). Hence the motion of any wave packet is also periodic with 
that same period. But there are important features of the general dynamics beyond this 
simple periodicity. To understand these it is instructive to study the motion of an initially 
localized Gaussian wave packet: 



V>(* = o) = £/*(o)|*> 



(26) 



with 



*(°)=(w) exp 



» 2 d 2 



4a 2 



+ ik 



(27) 



Fourier-transforming the time-evolution equation fe(t) = XV Kt,t'{t)ft'(0)i and utiliz- 
ing (|T0|), one obtains 



f(k,t) = Y / e- iMd K e/ ,(t)U(0) 

LP 



exp 



. A . eF dt\ ( eF dt 
i — , sin I — - — I cos kd + 



e 



-i(k+eF t)l' d 



(28) 



eF d "V 2 7 ~ V 2 
This equation is still exact. If we now assume that a/d^> 1, so that the spatial localization 
of the initial wave packet ( p^) is weak, then the summation over £' can approximately be 
replaced by integration, and we arrive at 



87TO" 
^2" 



2\ 1/4 



exp 



. A . eF dt\ eF dt 
i — ; sin I — - — I COS KC( + 



3 -o- 2 (fc-A:o+eFot) 2 



(29) 



eF d "' \ 2 J \ 2 
At this point we can make contact with two familiar notions. First, it can be seen that the 
center of the wave packet in fc-space is given by k(t) = ko — eF^t. Hence, 



jk{t) = -eF . 



(30) 



That is a special case of Bloch's famous "acceleration theorem" : the rate of change of the 
wave packet's center in fc-space is determined by the external force |26|]. Second, it should 
be realized that the explicit T^o^-periodicity of the discrete sum ( p8|) is apparently lost if 
one replaces summation by integration. But the wave vector k is defined, because of spatial 
periodicity, only modulo 27t/d, so this is also true of the argument (k — k + eF t) appearing 



11 



in the exponential of (|29|) , restoring the temporal periodicity with the Bloch period. The 
corresponding physics is that the wave packet is Bragg-reflected at the edges of the Brillouin 
zone. 

Now we are interested in the evolution of the wave packet in real space: 

h(t) = f[ +n/d dkf(k,t)e Md . (31) 

Z7T J-n/d 

Since we have already stipulated that a/d ^> 1, the momentum distribution f(k, t) is sharply 
localized around k(t). Therefore, it makes sense to expand the factor cos(kd + eFodt/2) in 
the exponential of around k(t). ^From a linear expansion one finds 

hit) « ' exp Uk - eF t)£d - i$(t) %)) 2 ) (32) 



with 



= ^Fl MM - eF dt) - sin(M)] (33) 



and 



£(t) = ^ [cos(k d - eF dt) - cos(A;orf)] . (34) 

Thus, the center of the wave packet, x(t) = £{t)d, oscillates in space with an amplitude that 
is inversely proportional to the dc field strength F , and with temporal period TbiocIi- This 
is the famous "Bloch oscillation" , that was first discussed by Zener |27] . 



A simple "semiclassical" explanation of these wave packet oscillations can be given as 
follows. First one introduces a continuous position variable x instead of the discrete lattice 
points Xe = id. From the "acceleration theorem" (|3C|) one knows k(t). Furthermore, the 
group velocity v(t) of the wave packet is determined by the dispersion relation ( f21~|) pertaining 
to the unperturbed lattice: 



= — sm(k(t)d) . (35) 

k(t) 



dk 

Assuming x(0) = and integrating, one finds for the position x(t) of the wave packet's 
center 
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x(t) = — — [cos(k d — eF dt) — cos(k d)] , (36) 
2eF 

which agrees exactly with £(t)d as derived above. The previous quantum mechanical cal- 
culation shows that this semiclassical reasoning is essentially correct as long as the spatial 
localization of the initial wave packet is weak, a/d ^> 1, so that the momentum distribution 
f(k, t) is sharply localized within a single Brillouin zone of width 27r/d. This assumption has 
not only been used when replacing sums by integrals, but also in the expansion employed 
to evaluate fl3~T|). It should be emphasized that the expansion around k(t) is not an expan- 
sion around a stationary point. It can, therefore, be justified only if the Gaussian in ( |29| ) 
is sharply peaked around k(t). That is the reason why we did not employ a second-order 
expansion, which would have resulted in a periodic oscillation of the wave packet width. 
But if the Gaussian is only moderately peaked, so that the second-order terms might be 
significant, the expansion around a non-stationary point is already questionable. 

However, the evolution of a wave packet that is sharply localized in space can easily be 
studied numerically. Fig. 1 shows an example for a = d, ko = 0, and eF$d = A/10. At 
time t = the wave packet is localized around the site £ = 0. The occupation probabilities 
|/f| 2 of the lattice sites, indicated by the open circles, have been connected by lines to 
guide the eye. The short dashes correspond to t = 0, the long dashes to t — TbiockI^-, 
and the full line represents the wave packet at t — Tsioch/2. After half a Bloch period, the 
maximum has moved by 9 sites to the left, whereas it would be 10 sites under "semiclassical" 
conditions. But even now the basic features of a Bloch oscillation are visible. The occupation 
probabilities at t — 3T^ oc / l /4 are the same as those for t = Tsioch/^, and at t = Tsioch the 
wave packet regains its initial shape. 

Fig. 2 depicts the same scenario for o = 5d; all other parameters are unchanged. This 
is a Bloch oscillation in the classical sense: the wave packet simply oscillates in space, with 
its shape remaining almost invariant. 

If the particle is initially localized entirely at a single site £ , i.e., if fi(0) = 5^, the time 
evolution is given directly by the propagator (p5|): 
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l/i(*)l s 



^(^sin(ew)) 



(37) 



In this extreme limiting case the wave function's center remains unchanged in time, but its 
width oscillates. 

IV. THE EFFECT OF COMBINED STATIC AND OSCILLATING FIELDS: 

DYNAMIC LOCALIZATION 

In order to solve the time-dependent Schrodinger equation idtip{t) = [Hq + Hi nt (t)]ijj(t) 
for the model described by (|IBD and (0), in the general case where H int (t) contains both a 
dc field F and an additional ac field Fx, it is helpful to introduce the vector potential 

F^ 

A(t) = -F t sin(wt) (38) 

and employ the gauge transformation 

#(t) = e' ieA(t)x ^j{t) , (39) 

where x denotes the position operator, x = ^i\i)id{l\. This restores spatial transla- 
tional symmetry: the Schrodinger equation for the new wave function ip(t) reads id t ip(t) = 
H(t)ip(t), where 

H(t) = (e- ieA{t)d \£ + 1){£\ + \£)(£+ l\e ieAit)d ) . (40) 

t 

Then the wave vector is a good quantum number in this gauge, even in the presence of the 
electric field. Observing that the Bloch waves \k (see (p0|)) obey the equation 

H(t) Xk =E(q k {t))xk (41) 

with 

q k (t) = k - eA(t) , (42) 

one immediately obtains solutions 
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V>fc(*) = El^> ex p(-^-* J o ^ E (lk(r))^ ■ (43) 
Inverting the transformation ([39|), one finds the wave functions 

Mt) = E \ £ ) exp(-iq k (t)£d - i J*dr E(q k (r))^ . (44) 

These are precisely the Houston states 123 for the driven tight-binding lattice. 

How are these states related to the Floquet states of the problem? After all, the Hamil- 
tonian Hq + Hi nt (t) is T-periodic, with T = 2tt/uj, so it might appear natural to search 
for T-periodic solutions u{t) to the eigenvalue equation [H + H int (t) — id t ]u{t) = eu{t). 
However, it must be recognized that the problem now actually contains two frequencies, 
the ac frequency u and the Bloch frequency u B [ och = eF d. As a consequence, the trans- 
formed Hamiltonian H(t) (see fl40|) ) is, in general, not T-periodic. Since H{t) contains both 
frequencies on an equal footing, one can then employ many-mode Floquet theory [p9|| . If 
^Bioch/u is a rational number, uo B i och /uj — p/q with p, q relatively prime integers, then H(t) 
is periodic in time with period T — q- 2h/uj = p ■ 2ix / uu B i och \ otherwise H(t) is quasiperiodic. 
The Bessel function expansion ( |i~0D readily yields 

E(q k (r)) = -- E Jr[^-)cos{kd+[eF d + ru}r) . (45) 

2 r^co \ U J 

If there is no integer n such that eF$d = nu, then the time integral over E(qk{r)) is an 
oscillating function with frequencies that are either multiples of the ac frequency u or of 
the Bloch frequency. In this case the wave functions (^) are two-mode Floquet states 
to the Hamiltonian ([4~0t). But this purely formal discussion is of limited physical value, 
since all quasienergies coincide (they are all zero, up to integer multiples of the respective 
frequencies), and do not depend on the field parameters. 

The situation is entirely different if usioch/^ — n for some integer n, — i.e., if 

eF d = nu , (46) 

so that the ac field is tuned to an "n-photon-resonance" with the Wannier-Stark ladder. 
Then H(t) is T-periodic, as is H + H int (t), so there are T-periodic Floquet states. Now 
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the expansion ( fl5f) contains the zero frequency term r = —n, which gives a linearly growing 
contribution to the total phase after time integration. But all other terms in the expansion 
of Jjdr E(qk(r)) are T-periodic. Hence, if ((46J) is satisfied, the quasienergies of the Floquet 
states (f|3|) are determined by the average growth rates (the zero frequency parts of (|45|)) of 
the phases [[3(j.fn|] : 

l r T 



e(k) = - f dtE(q k (t)) mod uj 
1 Jo 



TJ 

= (-l) n J n E(k) mod uj . (47) 

The gauge transformation fl39|) is a T-periodic operation only if ([ff]) is satisfied. Therefore, 
the Houston states ( pH| ) coincide with the Floquet states to H + H int (t) only in this resonant 
case, ipk{t) = Uk(t) exp(—ie(k)t). The quasimomentum k then remains a good quantum num- 
ber for the eigenfunctions «&(£), despite the presence of the dc field, so that ( f47D should be 
regarded as a quasienergy-quasimomentum dispersion relation. The particular significance 



of the resonance condition ([4g) has been discussed by Zak with emphasis on commuting 
time translations and electric translations. In the following, we will focus entirely on the 
resonant case. 

Whereas the Wannier-Stark energy eigenstates ( p2[) pertaining to dc fields are local- 
ized, an additional resonant ac field leads to extended quasienergy eigenstates that form 
quasienergy bands of finite width. Therfore, a resonant ac field can destroy Wannier-Stark 
localization, even if its amplitude is infinitesimally small. An arbitrary, initially localized 
wave packet ip(t = 0) can be represented as a superposition of Floquet states with a cer- 
tain momentum distribution f(k,0), and the time evolution of each component is simply 
determined by its quasienergy phase factor: 

f{k,T) = f{k,0)exp{-ie{k)T). (48) 

Assuming that the initial momentum distribution is sharply centered around k , the time- 
averaged group velocity of the wave packet is 



de(Jfe) 



dk 



-l) n ^jJ—)sm(kod). (49) 



2 "\ uj 

ko \ 
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Generally, therefore, the wave packet moves over the whole lattice, simultaneously spreading. 
Assuming the initial wave packet to be a Gaussian (|27|) with width <jq, the width er(sT) after 
s cycles is given, within the Gaussian approximation to the integrals defining the spatial 
Fourier transformation, by 

T2\ V2 



a{sT) = a 1 + 



Ad 2 T feF 1 d s 
4 



sT—^Jj cos(M) 



(50) 



Remarkably, the spreading of the wave packet is substantially hindered if k = ±n/2d, so 
that cos(kod) vanishes. The remaining slow spreading of the packet then is entirely related 
to terms beyond the Gaussian approximation. 

In the special case when eF^djio coincides with a zero j n ,s of the Bessel function J n , 
all quasienergies e(k) coincide. In that case (^(t)] 2 is strictly periodic in time with period 
T = 2tc/uj: on the average, the wave packet neither moves nor spreads. This has been called 



"dynamic localization" [p2| -|36[|. In view of the "collapse" of the quasienergy band |37| at 
the zeros of J n , dynamic localization can be interpreted as prohibited dephasing: according 
to P8|), all components of the wave function acquire exactly the same phase factor during 
one period T = 2tt/uj, so that the wave packet reassembles itself after every period, apart 
from an unimportant overall phase factor pj|| . 

There is a close correspondence between the approximate formula fllTD for the quasi- 
energies of the two-level system and its counterpart ([|7]) for the infinite lattice: if the centers 
of the two wells are separated by a distance d, then the dipole /i is approximately ed/2, so 



that the argument of the Bessel function in (17) becomes eFid/ui, as for the infinite lattice. 



But whereas the approximate formula (17) for the two-level system, or the formula (13) 



for the "renormalization" of the atomic ^-factor, is a perturbative result, the corresponding 
formula (f47D for the quasienergy band is exact. This is due to exact translational symmetry of 
the infinite lattice. It is quite remarkable that the modification of the energy level structure 
by a Bessel function, which appears quite naturally in the context of infinite lattices, still 
remains a good approximation for finite lattices, and even survives in the limiting case of a 
"lattice" with only two sites [38]. 
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If one formulates the problem of the interaction of particles in spatially periodic potentials 
with temporally periodic fields in terms of Floquet states and quasienergy bands, one is 
adopting a point of view quite similar to the "dressed atom" -picture PJItj in atomic physics: 
the ac field is not regarded as a perturbation. Rather, the lattice and the ac field together 
form a new object, the "dressed lattice", and this object has its own quasienergy bands. 
Exactly as an ac field can "renormalize" an atomic g-factor, it can also "renormalize" the 
energy band structure of a lattice. In the following sections we will show that this point of 
view is particularly useful for understanding effects caused by lattice imperfections in the 
presence of strong ac fields. 

V. ISOLATED DEFECTS: MANIPULATING LOCALIZATION LENGTHS BY AC 

FIELDS 

Now we introduce an isolated defect into the lattice and study the dynamics governed 
by the Hamiltonian, H(t) = H Q + H int (t) + H^efect, where the defect Hamiltonian 



changes the energy of the atomic state |0) by Uq. H and H int (t) are again given by fllBp 



the behavior of the Floquet state localized around the defect. 

If we were to consider the simpler problem of a time-independent perturbation V of the 
lattice Hamiltonian Hq, we could employ the usual resolvent operator formalism: the energy 
eigenvalues of the perturbed Hamiltonian H = Hq + V are given by the poles of the resolvent 



H defect = vo\Q)(Q\ 



(51) 



and (|I5D, respectively. We assume that the ac field is resonant, eF d = nu, and investigate 




(52) 



in the complex z-plane, where 




(53) 
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is the resolvent of the unperturbed system. In a representation of spatially localized states 
the resolvent equation (|5^) becomes the Dyson equation 

Rjt(z) = R ( $(z) + ]T Rf}{z)V re Rei{z) . (54) 
i'i' 

If the perturbation V has only a finite number of nonvanishing matrix elements Vyf, then 
the right-hand side of ( |54] ) involves Re>e(z) only for a finite number of sites £'. Therefore, af- 
ter ( Ell ) has been inverted for that subset of resolvent operator matrix elements, the equation 
itself gives Rje{z) for all sites j. 

The task now is to extend this formalism to the case where there is an additional time- 
periodic external field [39]. To this end we regard the lattice together with the driving field 



as the unperturbed system; the perturbation is just the defect potential. In order to solve 
the eigenvalue equation (H), we have to compute the resolvent of the quasienergy operator 

H = H(t) - id t (55) 

in the extended Hilbert space of T-periodic functions. 

A complete set of spatially localized states in that space is given by the functions 

\£, m) = \£) exp(ieA(t)£d + imtut) (56) 

with the vector potential A(t) (see (^)). With respect to this set we define Green's functions 

G je (n,m) = {{j,n\ ^ \£,m)) , (57) 

where the double brackets indicate the scalar product (f|). We no longer explicitly indicate 
the dependence of these Green's functions on z. The Dyson equation becomes 

Gji(n, m) = G$(n, m) + ]T G$(n, n') V fe G t , e (ri, m) , (58) 

j'i'n' 

where Gjf(n,m) are the unperturbed Green's functions pertaining to the operator (H + 
Hintit) — id t ), which does not include the time-independent defect potential V. Compared 
to (|5"4"D, the summation now also includes the additional Fourier index n' . 
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In order to compute the unperturbed Green's functions, we employ the spatially extended 
quasienergy eigenfunctions (with integer s) 

\k,s} = u k (t)e^ t 

= 1-0 exp (-iq k (t)id - i J*dr [E{q k {r)) - e{k)} + isutj (59) 



constructed from the Houston states fl44|). Curly brackets are used to distinguish these 
eigenfunctions from the localized functions fl56|). We then have 



€${n,m) = £<0>IM}} — {{k,s\£,m)) . (60) 

M z-e[k)-suj 

The matrix elements can now be expressed as 

((j,n\k,s}} = e-^ d F s _ n , (61) 

where 

Fs ' n = T £ dtel(S ~ n)uJtex V / * dT - e(A:)]) . (62) 

Obviously, the formalism considered here can be of practical use only if it is possible to 
neglect all but a finite, small number of Fourier components. Such an approximation is 
justified in the case of high frequencies: the expansion (cf. (ft5|)) 

ft , r „. ... A ^ T ( eFid\ sin(kd + (n + r)ut) — sin(kd) , . 
/„ <W»M> "«(*)] = - Y ^ J r (— J L_1_J LJ (63 ) 

shows that the integral on the left hand side is — at most — of order A/a;; the oscillatory 
character of the function E(k) makes it, in general, much smaller than this upper limit, as 
shown by the explicit expression on the right hand side. Therefore, for frequencies u larger 
than the bandwidth A a reasonable approximation to ( |62]) is 

F s - n ^5 s , n , (64) 

and the unperturbed Green's functions become 

,(0), , d r/ d „ (pw-fi 



Gy(n,m)tt—l dk — 5 nm . (65) 

2n J-Tt/d z - e(k) -nuo 
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In this high-frequency limit the Green's functions are diagonal in the Fourier index. In the 
following we will restrict ourselves to this case, and suppress the redundant Fourier index 
altogether. 

After inserting the quasienergy-quasimomentum dispersion relation 



where 



e(k) = (-l)("+ 1 )jJ^j^cos(A ; d) = ycos(fcrf), 



(66) 



(67) 



denotes the width of the quasienergy band (up to a sign), the unperturbed Green's functions 
can be calculated with the help of the residue theorem: 



G K 



(0) 



4 C. 



l^-j'l+i 
o 



ji Wl- Co 2 



2 ' 



where 



_2z Uz 2 



The sign has to be chosen such that £o falls inside the unit circle. 
For an isolated defect of the type given by fl5T]) we now have 



(68) 



(69) 



(70) 



and the Dyson equation (|58 ) becomes in the high-frequency limit 



(71) 



Specifying j = 0, solving for Goe, and reinserting, one obtains 

Uif — u 



1 - ^ ' 



r 00 



The pole which lies inside the unit circle, |£o| < 1, is found to be 

2u 



Co = sign(z/ /W) 



W 2 



+ 1 - 



w 



(72) 



(73) 
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Hence, the quasienergy of the state supported by the defect is given by 



W ( 1 \ W 2 
Sdefect = -j (Co + -y J = sign(z/ )yz/2 + — . (74) 

Thus, while the width W of the quasienergy band oscillates when the parameter eFid/u is 
varied, the quasienergy of the defect state follows these oscillations, maintaining a positive 
distance from the band. 

The oscillations of the band width have a pronounced effect on the spatial density of 
the defect state. In the high-frequency regime, the probability p(£) for the defect state to 
occupy the £-th site is determined by the residue of the diagonal Green's function Gu in the 



complex quasienergy plane p9fl . An elementary calculation yields 

2\e\ 

(75) 



v 2 + W 2 /A \V W 2 



w 



Thus, the occupation probability falls off exponentially from the site I = where the defect 
is located. The inverse exponential decay length, measured in multiples of the lattice spacing 
d, is given by 



L\ _1 / /4z/ n 2 2z/ 

-21n \/— ^ + 1- 



W 



(76) 



d) \\w 2 

Obviously, L~ l increases monotonically with Iz/q/H^I, the ratio of the size of the perturbation 
and the quasienergy band width. Therefore, it is possible to manipulate the spatial extension 
of defect states by tuning the parameters of the driving field. When eFid/u approaches a 
zero of the Bessel function J n , while eF d = nu, the defect state should be concentrated 
wholly at the anomalous site This possibility of controlling localization lengths by 

external fields is one of the most interesting predictions of the "dressed lattice" -picture. 

To check the validity of our conclusions even away from the high frequency limit, we 
perform numerical computations for a finite lattice with iV = 101 sites, with a defect placed 
in the center. The wave functions vanish at the chain ends. Fig. 3 shows the quasienergy 
spectrum for u/A = 1.0 and eF^d = u, i.e., for a "one-photon resonance". The defect 
strength is z/ /A = 0.1. The quasienergy band structure is determined by the Bessel function 

22 



Ji, the "collapse" of the band at = 3.83171 and ji^ = 7.01559 is indicated by the arrows. 
As expected, the quasienergy of the defect state on top of the band follows the oscillations 
of the band width. 

To investigate the degree of localization of the defect state, we compute its quasienergy 
eigenfunction 

N 

u defect (t) = Y,c ( i defect \t)\£) (77) 
and calculate the time- dependent inverse participation ratio 

N 

P, e/erf (t)=El4 de/eCt) W| 4 - (78) 
1=1 

If the defect state were uniformly extended over all N sites, Pdefect{t) would vanish as 1/N, 
whereas it approaches unity for a completely localized state. Fig. 4 shows numerical results 
for the same situation as considered in Fig. 3. The values of the scaled ac field strength 
eFid/u are (short dashes), 5.0 (long dashes), and (full line). The density of the 
defect state exhibits a periodic oscillation that leads to the observed time-dependene of the 
inverse participation ratio. These oscillations become negligibly small for higher frequencies; 
Pdefect then becomes independent of time. But clearly, even for u — A the localization of 
the defect state is much more pronounced at the zeros of the Bessel function than at other 
parameters. 

In order to test the prediction of ([75]) directly against the numerical data, we compute 
the participation ratio as a function of eF\d/uj at t — 3T/4, when the ac field vanishes. The 
result is shown as the full line in Fig. 5. The broken line, on the other hand, results from 
the evaluation of the approximate formula (|75|) : Pdefect ~ X^G ^)) 2 - The two lines differ for 
low field strengths, as the approximate expression goes to unity when F\ vanishes, whereas 
the exact result approaches a smaller limit, representing a state spread over more than the 
single impurity site. However, the agreement is strikingly good for large ac amplitudes. We 
conclude that the high-frequency approximation can give good results even if the frequency 
uj is not really large compared to the bandwidth A. 
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It should be stressed that in the high-frequency regime uo > A where the approxima- 
tion fl64"|) is valid the effect of the driving field is fully described by the quasienergy spectrum: 
the driven lattice with quasienergy band width W becomes equivalent to an undriven lat- 
tice with that same energy band width. This observation stets the stage for the following 
discussion of Anderson-type localization in the presence of ac fields. 



VI. RANDOM DISORDER: CONTROLLING ANDERSON LOCALIZATION 

Having dealt with the effects of a single defect, we now consider randomly disordered 
lattices: 

H(t) = H + Hi nt (t) + H random (79) 



where H and H int {t) are given by flT8|) and (|19"D , as before, and H ran dom introduces site- 
diagonal disorder: 

H randvm = '£v e \£)(£\ . (80) 

i 

For specificity we take the probability distribution p[y) for the random energies U£ as 

p{u) = — . 1 (81) 

ITU / t \2 

" "max / v \ 

V \ Umax J 

for \u\ < Umax, and zero otherwise, though the precise choice of p(u) does not greatly 
affect the conclusions below. Without external fields, the time independent Hamiltonian 
Hq + H ran do m defines a paradigmatic system for the study of Anderson localization |40|-|43|. 
It admits only localized states, for arbitrary positive disorder strength v max . The typical 
localization length L is given by the square of the ratio of the bandwidth A and u v 

A ^ 2 



d . (82) 



Thus, if the disorder strength becomes comparable to the unperturbed bandwidth, the 
extension of typical energy eigenstates is only of the order of a few sites. 
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Consequently, the fact that in the presence of resonant high-frequency ac fields the 
amplitude-dependent quasienergy band width W takes over the role of the energy band 
width A leads to an interesting prediction |45|]: It must be possible to manipulate the 
degree of Anderson localization by tuning the dimensionless ac field strength z = eFid/uu, 
provided the ac field is resonant with the Wannier-Stark ladder, i.e., eF^d = nuo for some 
integer n. To demonstrate this effect, we compute the Floquet eigenf unctions for a finite 
disordered lattice with N = 101 sites, 

N 

u m (t) = J2c ( r\t)\£), (83) 
1=1 

and then determine the average inverse participation ratio P: 

1 N 

f^E |c^(3T/4)| 4 . (84) 

£,m=l 

As in the previous case of an isolated defect, the choice of the particular moment t = 3T/4 
is motivated by the fact that the ac field vanishes at that time (see (|T9|)), but it is without 
special significance. In the high-frequency limit, where the time-dependence of the Floquet 
states becomes weak, P is almost independent of the particular time at which it is evaluated. 

Fig. 6 shows P versus z = eFid/u for n = 1 and u — A, with v max //S. = 0.01, 0.05, 
0.10, and 0.20. The effect of the collapse of the quasienergy band is quite dramatic: when z 
approaches a zero of J\, the Floquet states are localized almost entirely at individual sites, 
as testified by values of P close to unity. 

We thus have to distinguish clearly the following effects: A dc field applied to a tight- 
binding lattice results in Wannier-Stark localization; every wave packet can be represented 
as a linear combination of localized energy eigenstates. An additional ac field, tuned exactly 
to an "n-photon resonance" with the Wannier-Stark ladder, leads to extended quasienergy 
eigenstates, and an initially localized wave packet will, in general, move throughout the 
lattice and spread, unless eFid/u equals a zero of J n . At these special parameters there 
is no quasienergy-quasimomentum dispersion, and the wave packet, although built up from 
extended states, remains "dynamically" localized. If there is additional random disorder, 
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the quasienergy eigenstates exhibit Anderson localization, with amplitude-dependent local- 
ization lengths. Anderson localization is most pronounced at the very same parameters 
where there would be dynamic localization in an ideal lattice. To illustrate the profound 
difference between dynamic localization and field-induced Anderson localization, we show in 
Fig. 7 the quasienergy spectrum corresponding to the case of strongest disorder considered 
in Fig. 6, v max j A = 0.20. At the zeros of J\, where the quasienergy band of the ideal lattice 
would collapse — a necessary condition for dynamic localization [RD| - the quasienergies 
of the disordered lattice are scattered over a range of 2v max . The components of an initially 
localized wave packet will, therefore, dephase in the course of time, in contrast to the ideal 
case. Nevertheless, the wave packet will remain localized, since all the Floquet states are 
individually localized themselves. 

In short: a dc field gives rise to localized energy eigenstates, an additional resonant ac 
field then leads to extended quasienergy eigenstates, and random disorder results in An- 
derson localization of these quasienergy states. Thus, localization in randomly disordered, 
resonantly driven tight-binding lattices has nothing to do with the Wannier-Stark localiza- 
tion that would result from a dc field applied to an ideal lattice. 

To substantiate this statement, we confirm the scaling behavior expected for genuine 
Anderson-type localization. The relevant dimensionless parameter of the problem should 
be the ratio of the quasienergy band width and the disorder strength. If this hypothesis is 
correct, then lattices with different disorder strengths must yield the same average inverse 
participation ratio P, if their values of AJ n (eFid/u) / u max coincide. In Figs. 8, 9, and 10 we 
plot P for n = 1, n — 2, and n = 3, respectively, again for lattices with N = 101 sites and 
u/A — 1.0. The disorder strength varies from v max j A = 0.01 to 0.10, in steps of 0.01. As 
can be seen, a certain characteristic ac field strength is necessary to effectively counteract 
the localization found for vanishing ac amplitude. As a measure of this characteristic field 
strength we take that value z\/i of the dimensionless amplitude z = eFid/uj for which P(z) 
is reduced to half the "Wannier-Stark" -value: P(z\/2) = P(0)/2. Since P should depend 
only on the ratio Aj n (z)/v max , there should be a relation [ffl 
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AJ n (z 1/2 ) = c n v max , (85) 

with numbers c n of order unity. The n-dependence of these numbers results from the fact 
that the degree of localization for z = depends on n, i.e., on the strength of the dc field. 
As long as the characteristic ac amplitudes remain weak, zy 2 < n/2, one can approximate 
the Bessel functions J n by their lowest order terms, J n (z) ~ (z/2) n /n\, to get 

v ~~ \ l / n 



z 1/2 « 2 \ n\ c n —£-J ■ (86) 

Provided the disorder is so weak that the characteristic ac amplitudes also remain suffi- 
ciently small, these amplitudes should scale approximately with the n-th root of the disorder 



strength j46|| . 



Fig. 11 confirms these hypotheses. The boxes denote the characteristic amplitudes for 
n — 1, 2, and 3, as determined from the preceding three figures. The dashed lines follow 



from the approximate relation (86), with c\ = 1.074, c 2 = 1.183, and C3 = 1.189. The 
agreement with the numerical data points is very good for n — 1 and n — 2, though there 
are clear deviations for n = 3. But these deviations merely indicate that the lowest-order 
approximation to J3 is already insufficient: if the full Anderson relation (|85| ) is evaluated 
without that approximation one obtains the solid line, which gives a truly excellent fit to the 
numerical data. This impressive agreement between the numerical results and the Anderson 
relation ( jSq) confirms that we are not dealing with remnants of Wannier-Stark localization 
induced by the dc field, but rather with Anderson localization resulting solely from disorder. 

^From a dynamical point of view, it is of interest to study directly the effect that Ander- 
son localization of quasienergy states has on the time evolution of an initially localized wave 
packet. To this end, we propagate a wave packet given at t = by a Gaussian (|27D in a lattice 
with N = 1001 sites. The field parameters are u/A = 1.0, eF d/u = 1.0, and eFid/u = 0.3. 
The diagonal disorder is again distributed according to (|8"ID , with u max /A = 0.1. The initial 
width of the packet is a = 5d, and the initial momentum k = it / (Ad) . Fig. 12 depicts the 
deviation of the wave packet's center, x = Y^£^d\fi(t)\ 2 , from the initial position x Q = 0, 
whereas Fig. 13 shows the evolution of its width a, both for a period of 1000 cycles T = 2n /oj 
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of the external field. As expected, the packet does not move on the average, the typical ex- 
cursion length being merely one or two lattice constants, and its width does not increase 
appreciably. If there were no disorder, the packet would move by more than 300 lattice sites 
to the left within 1000 ac cycles, and its width would grow as indicated by the dashed line 
in Fig. 13. 

As in the case of an isolated defect, the identification of a periodically driven lattice 
of bandwidth A with an undriven lattice of "renormalized" bandwidth AJ n (eFid/u) is 
essentially correct in the high-frequency regime, to/ A > 1. Nevertheless, localization effects 
can efficiently be manipulated by ac fields even if this condition is not satisfied, provided the 
dc field alone is strong enough to cause substantial Wannier-Stark localization. Fig. 14 shows 
an example for toj A = 0.1 and ei<o<i = 10 u, i.e., for a "10-photon resonance". The lattice has 
101 sites, and the disorder strength is v max j A = 0.1. The average inverse participation ratio 
P remains almost constant for weak amplitudes, then shows a remarkably sharp crossover 
from "strong" to "weak" localization, and exceeds 0.5 again in the vicinity of jio,i ~ 14.5. 
At the second positive zero of Jio, i.e., at jio,2 ~ 18.4, localization is less pronounced. 

Fig. 15 shows the quasienergy spectrum for this case. Since u <C A, the quasienergy 
band of the ideal lattice does not fit into a single quasienergy Brillouin zone, but overlaps 
itself. Because of the loss of translational invariance due to disorder, this leads to a dense 
mesh of avoided crossings when the quasienergies of the random lattice are plotted, e.g., 
versus the ac amplitude. But even though the quasienergy spectrum appears quite intricate 
in the low-frequency case, there remains a clearly visible connection between the behavior 
of the eigenvalues and that of the inverse participation ratio. 



VII. INTERBAND TRANSITIONS 

The discussion so far has completely neglected transitions between different energy bands. 
But even if there is only a dc field, interband transitions can give rise to important and 
interesting changes of the dynamics ||2T|J4"T| , |4*8|| . In order to study the role of interband 
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transitions in the presence of combined ac and dc fields, we now investigate the two-band 



model defined by the Hamiltonian [25 



H{t) = H^it) + H {2 \t) + H mterband (t) . (87) 

For either j = 1 or j = 2 the Hamiltonian H^'(t) describes a particle in a driven tight- 
binding lattice with a single band of width Aj: 

H®® = (-iyB. y: \£,j)(£,j\ + eF(t)dY: \tj)t(£,j\ 

+(-i) i ^E(l^ 1 .Jl( < .jl + l^'>^ 1 .il) • ( 88 ) 

F(t) is the total external field, 

F(t) = F + F 1 cos(ut) . (89) 

If both Fq and F\ vanish, then D denotes the energy difference between the band centers. 
The interaction between the two bands is modelled by 

H mt erban d (t) = eF(t)X 1>2 £ (\£, l)(£,2\ + \£,2)(£, 1|) . (90) 

e 

Here is the matrix element of position x between the Wannier states of the two bands 
centered on a single site. If we approximate it by the dipole matrix element between the 
lowest two energy eigenstates of a particle in a box of width d, as in j47|, we have = 



— 16d/(97i 2 ). We neglect, as usual, dipole matrix elements between Wannier states centered 
on different sites. 

If the ac field is resonant with each of the two interpenetrating Wannier-Stark ladders 
originating from the two bands [25], eF d = nu, then there are quasienergy bands with 
finite widths. For vanishing interband interaction, = 0, the normalized quasienergy 
eigenf unctions for a lattice with N sites are (neglecting finite size effects) 

ufit) = "4 E %3) {-iq k (t)£d - i /W [E^(q k (r)) - e®(k)]) , (91) 
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as in (13), where 



EV\k) = (-I)'- + (-iy — cos(kd) (92) 

and 

e®{k) = (-l)'f + (-l)°' +n) ^-4(^) «>s(faQ (93) 

are the energy and quasienergy eigenvalues, respectively. 

For nonzero X 12 , and resonant ac fields, the two-band model (|87|) is equivalent to a 
system of two-level systems labelled by the wave vector k. If (Ai + A 2 )/2 < cj, which will 
be fulfilled in many situations of interest, a high-frequency approximation analogous to fl64|) 
yields the matrix elements (r, s integer) 

(94) 

The discussion of section 2 then shows that, besides the resonance eF d = nuo required for 
the emergence of quasienergy bands, there is a second, less obvious resonance: if also 

2eF Xi i2 = mu (95) 

for some integer m, then the quasienergy bands will be most strongly affected by H interband {t). 

One can also arrive at this second resonance condition (|95|) by a different argument. If 
both Ai = and A 2 = 0, then (|87|) reduces to a system of uncoupled two-level systems 
Hfo(t) that are labelled by the site index t: 

Hf s (t) = ja z + eF(t)£d 1 + eF(t)X^ 2 a x . (96) 

All of these are unitarily equivalent to H^(t), which is just ([TJ|) with A and /1 replaced by D 
and eXi t 2, respectively. Thus, if (|95|) holds, the quasienergies are determined approximately 
by ([TJ]) for strong fields or high frequencies, for each I. For nonvanishing Aj this ^-degeneracy 
is removed, and the quasienergy bands obtain finite widths. 

We illustrate these qualitative considerations by numerical computations. First, Fig. 16 
shows the quasienergies for the two-level systems (Off) with eF Q d = to and Xl i2 = —16<i/ (97r 2 ), 
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so that (p5|) is not satisfied. The energy level separation is D = 5.5 id, hence uj/D < 1. The 
behavior of the quasienergies in this low-frequency case, including the appearance of the 
sharp avoided crossing at eFid/u pa 17, can be understood with the help of the arguments 
given in ref . [|30| . 



Fig. 17 then depicts the quasienergy spectrum for the two-band model fl57|) with the 
same parameters D/u, eFod/u, and X\^/d. The number of sites is AT = 20; the bandwidths 
are Ai = 1.2 to and A 2 = id. The centers of the two quasienergy bands coincide precisely 
with the quasienergies of the associated two-level system, cf. Fig. 16. Since we have a 
"one-photon-resonance", the oscillations of the bandwidths are determined by J\{eF\d/ui). 
Arrows on top of the figure indicate values of the scaled ac amplitude eF\d/uj that coincide 
with a zero of J\. At these zeros, both bands collapse. In addition, there is strong level 
repulsion at the borders (e/u = ±1/2) and in the center (e/u = 0) of the quasienergy 
Brillouin zone. 

When the interband separation D is decreased, the dynamics become more complicated. 
Fig. 18 shows quasienergies for D = 1.6 uj; all other parameters are as in Fig. 17. Some 
of the band collapses now become imperfect, similar to the observations in ref. j|7|], and 
additional band narrowings occur. 

In order to satisfy the second resonance condition (|95p together with the basic condition 
eFod = noo, it is necessary that 2X\^jd = m/n be a ratio of two integers. Of course, it 
is always possible to find rational approximants m/n arbitrarily close to any given value 
of 2X^2/ d, but the most relevant cases are those where both n and m are small. Fig. 19 
depicts the quasienergy bands for eF\d = u and 2X\^/d = —1. The other parameters are 
D = 3.5 U, Ax = 1.2 id, and A 2 = id. Again, arrows on top of the figure indicate zeros of 
J\. As expected, the quasienergy spectrum now appears quite different from those shown in 
Figs. 17 and 18. In particular, there are no band collapses at high ac amplitudes. 

The first avoided quasienergy band crossing that occurs in Fig. 19 in the center of the 
Brillouin zone is shown enlarged in Fig. 20. Such avoided crossings are of great significance: 
at an avoided crossing, Fourier spectra of Floquet functions peak strongly around a certain 
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harmonic of the ac frequency [j30| , [l9|| . This observation has led to the suggestion that a 
semiconductor superlattice, irradiated by a strong far-infrared laser field under conditions 
which correspond to an avoided crossing of its quasienergy minibands, might emit radiation 
preferentially at some odd integer multiple of the driving frequency uo [ f49f . 

The merit of the introduction of quasienergy bands is an enormous conceptual simpli- 
fication. If one uses the energy eigenstates, i.e., the unperturbed Bloch waves, to describe 
the dynamics of an electron in a lattice under simultaneous ac and dc bias, one first has to 



account for Zener tunneling that occurs in static fields ||27|| . Next, the oscillating field leads 
to Rabi-type oscillations between the energy eigenstates. The systematic investigation of 
the effects that additional lattice disorder would have then appears almost hopeless. The 
Floquet picture, on the other hand, is much more economical, since it fully incorporates 
both external fields in the basis of quasienergy eigenfunctions. The influence of disorder 
can then be discussed with the help of the quasienergy spectrum: the degree of localization 



depends on the quasienergy band widths [50 



VIII. DISCUSSION AND OUTLOOK 

Let us summarize: The quantum mechanics of particles in periodic potentials that in- 
teract simultaneously with static and with resonant, oscillating electric fields can most con- 
veniently be studied by introducing quasienergy states instead of the usual Bloch waves. 
The widths of the quasienergy bands depend strongly on the field parameters, and can even 
approach zero. This band width renormalization bears close mathematical and conceptual 
similarity to the well-studied renormalization of atomic g-factors |]-§|]- If there is lattice 
disorder, the ratio of disorder strength and quasienergy band width determines the localiza- 
tion lengths of quasienergy eigenstates, exactly as the ratio of disorder strength and energy 
band width determines the localization of energy eigenstates in the static case. By adjusting 
the ac amplitude, it is possible to switch from the regime of "weak" localization, where the 
localization lengths are well beyond the length of the entire lattice, to the regime of "strong" 
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localization, where the quasienergy states are essentially localized at individual sites. 

We have restricted the discussion to model systems with nearest neighbor interactions 



only, but it is not difficult to drop this restriction [ST|. For example, consider the case where 
there are also next nearest neighbor couplings. Then the dispersion relation of the undriven 
lattice contains an additional term proportional to cos(2kd), so that quasienergy bands of 
finite width occur not only for eF^d = nut, but also for eF$d = nu/2. The appearance 
of additional bands reflects an n-photon-resonance between "atomic" states that are next 
nearest neighbors. 

It has already been mentioned that semiconductor superlattices are possible candidates 
for an experimental test of the ideas outlined in this paper. For such superlattices the nearest 
neighbor approximation is appropriate, and there is always some disorder in these artificially 
grown mesoscopic systems. It must be emphasized that an electron in a semiconductor expe- 
riences not only the external field, but also f52| , |53| a field from induced polarization. However, 
it can be shown theoretically that dynamic localization in ideal superlattices should persist 
even in the presence of Coulomb interactions [|54|,[55|], and signatures of dynamic localiza- 
tion have actually been observed in experiments performed at UCSB ]15|]. Therefore, we 
may assume that our models, although highly idealized, still capture a significant part of the 
physics of real superlattices. A recent analysis |55| of a realistic full three-dimensional model 
of an anisotropic semiconductor, including exciton formation as a consequence of coulomb 
interactions, demonstrates that the dynamic localization phenomenon is accompanied by a 
change in the effective dimensionality of the superlattice exciton. 

Entirely different examples of particles in periodically driven, one- dimensional periodic 
potentials are provided by ultracold atoms in modulated standing light waves. If atoms 
of mass M move in a monochromatic standing light wave with a frequeny that is detuned 
sufficiently far from any atomic resonance, and if the standing wave position is periodically 
modulated with an amplitude AL and angular frequency u m , then the effective Hamiltonian 
can be written as |56|-p9|] 
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H (t) = Si ~ % cos[2A; L (x - AL sw(u m t))] , (97) 

ZlVl o 

where Q denotes the effective Rabi frequency, and ki is the light wave number. This Hamil- 
tonian is unitarily equivalent, and quasienergetically isospectral, to 

^ {t) = 2M ~ I C0 < 2kLX ^ ~ Xx sin ^) - 4^2" ' ( 98 ) 



in 



with 

A = Muj 2 m AL . (99) 

The corresponding classical Hamiltonian system exhibits chaotic dynamics. The above 
model system has been used to study, both experimentally and theoretically, the quan- 
tum mechanical suppression of classical phase space diffusion |56|-|59j. (Incidentally, this 
suppression of classical diffusion has also been termed "dynamical localization" . It must be 
recognized, however, that it is in no way related to the "dynamic localization" discussed in 
section 4.) 

Clearly, ultracold atoms in standing light waves are higly attractive candidates for study- 
ing quantum dynamics in periodic potentials, because of the conceptual simplicity of the 
experimental set-up and the high degree of control over the parameters ||56|| . In fact, one 



can not only tailor the initial momentum distribution of the atoms at will, but also turn 
the periodic potential off easily, thus obtaining access to the momentum distribution of de 
Broglie waves that have evolved under the influence of that potential []60 |. Moreover, there 



is no dissipation or scattering from defects, and the interaction between the neutral atoms 
can be negelected. These unique features have been exploited in a very recent beautiful 



experiment [60]: it has been demonstrated that ultracold Cesium atoms prepared in the 



ground energy band of the potential created by a standing light wave perform Bloch oscil- 



lations when they are driven by a constant force. Ben Dahan et al. [60] were able to deduce 
the velocity of the oscillating atoms as function of time and obtained perfect agreement 
with the theoretical prediction, thus visualizing for the first time what Bloch and Zener had 
envisioned decades ago. 
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In the light of these achievements, other possibilities come to mind readily. If the current 
experiments with periodically forced atoms j|56| , |5"8| can be extended into the single-band 
tight-binding regime, it should become feasible to study dynamic localization (in the sense 
of section 4) of atomic de Broglie waves — i.e., spatial localization of ultracold atoms in 
standing light waves, caused by collapses of their quasienergy bands. Another interesting 
question is suggested by eq. (|50|): can one transport an atomic wave packet with suitably 
prepared initial momentum in the field of a standing light wave (almost) without dispersion? 
It appears quite likely that the already successful combination of atomic, optical, and solid 
state physics will unravel further important effects. 
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FIGURES 

FIG. 1. Time evolution of a sharply localized wave packet ( p7|) in a tight-binding lattice in the 
presence of a static electric field, eFod = A/10. The initial localization length a is a single lattice 
period: a = d. The site populations \fe(t)\ 2 are connected by lines to guide the eye. Short dashes: 
t = 0, long dashes: t = Tsioch/^, full line: t = T^iochl^^ where Tsioch = ^/(eF^d) denotes the 
Bloch period. 



FIG. 2. Same plots as Fig. 1, but for a wave packet not so sharply localized: a = 5d. 

FIG. 3. Quasienergy spectrum for a finite lattice with N = 101 sites, and an isolated defect 
in the center. The ac frequency is u>/A = 1.0, and eF$d = uj, corresponding to a "one-photon 
resonance". The predicted collapse of the quasienergy band at = 3.83171 and j\p = 7.01559 
is indicated by the arrows. The defect strength is z/q/A = 0.1. The quasienergy of the defect state 



appears above the quasienergy band, cf. (74). 



FIG. 4. Time-dependent inverse participation ratio Pdefect(t) (see (|7"8|)) for the same situation 
as considered in Fig. 3, for eF\d/uj = (short dashes), 5.0 (long dashes), and jip (full line). The 
origin of time in this figure, t = 0, corresponds to t = 3T/4 in the Hamiltonian (|l~9|), when the ac 
field vanishes. 

FIG. 5. Numerically determined inverse participation ratio Pd e / ec t(3T/4) (full line) versus 
scaled ac amplitude, for the same situation as considered in Fig. 3, compared to the high-frequency 



approximation (75) (dashed line) 



FIG. 6. Average inverse participation ratios ( |S4"|) for a driven, disordered lattice with N = 101 
sites, as functions of the dimensionless ac amplitude eFid/u, showing the remarkably sharp 
crossover from "weakly" to "strongly" localized states. The ac frequency is lo/A = 1.0, the dc 
field strength eF^d/oj = 1.0. The disorder distribution is given by (|Sl|), with u max /A = 0.01, 0.05, 
0.10, and 0.20. The weaker the disorder, the sharper the spikes. 
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FIG. 7. Quasienergy spectrum for a driven lattice with disorder strength u max /A = 0.20; all 
other parameters are as in Fig. 6. 

FIG. 8. Average inverse participation ratios for the same lattice as considered in Fig. 6, again 
for the "one-photon resonance" eF^d = 1 • oj. The disorder strengths vary from v max j A = 0.01 
(lowest curve) to 0.10. 



FIG. 9. Same as Fig. 8, for eF d = 2 -to. 



FIG. 10. Same as Fig. 8, for eF d = 3-uj. 

FIG. 11. Boxes: characteristic ac field strengths Z1/2 versus disorder strength for n = 1, 2, 
and 3, as determined from Figs. 8-10. Dashed lines: approximate relation (|86| ) with c\ = 1.074, 
C2 = 1.183, and C3 = 1.189. Full line: prediction of the full Anderson relation (|85|) for n = 3. 



FIG. 12. Time evolution of the center x of a wave packet initially given by a Gaussian (27) in a 
disordered lattice. At t = the packet is localized at the center of a lattice with 1001 sites. The ac 
frequency is oj/A = 1.0, the dc field strength eF^d/u) = 1.0, and the ac field strength eF\djuj = 0.3. 
The disorder strength is v max j A = 0.1. Without disorder, the packet would move to the left by 
328 sites within 1000 cycles T = 2tt/uj. 

FIG. 13. Time evolution of the width of the same wave packet as considered in Fig. 12 (full 
line). The initial width is a = 5d. The evolution of the width in the absence of disorder is shown 
by the dashed line; a would reach 33.2 d after 1000 cycles. 

FIG. 14. Average inverse participation ratio for a low-frequency field: uj/A = 0.1, with 
eF$d = 10 oj. The lattice has 101 sites; the disorder strength is v max j A = 0.1. The sharp spike is 
at the first zero of J\q, where eF\d/u) ~ 14.5. 



FIG. 15. Quasienergy spectrum corresponding to the situation shown in Fig. 14. 
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FIG. 16. Quasienergies for the two-level systems (^g) with D/uj = 5.5, and eF$d = u. The 
interband matrix element is X\^/d = — 16/(97r 2 ), so that the resonance condition (^) is not 
satisfied. These are also the quasienergies for the two-band model with Ai = A2 = 0. 

FIG. 17. Quasienergies for the two-band model ( |S7| ) with 20 sites. The parameters D/uj, 
eFod/uj, and Xi^/d are as in Fig. 16; the bandwidths are Ai/a> = 1.2 and A^/uj = 1.0. Arrows 
on top of the figure mark values of the scaled ac amplitude eF\d/uj equal to zeros of the Bessel 
function J\. 



FIG. 18. Same as Fig. 17, with reduced interband separation D/uj = 1.6. 

FIG. 19. Quasienergies for the two-band model (jSTj) with D/uj = 3.5, Ai/u = 1.2, A2/0J = 1.0, 
and eFod = u. Now the interband matrix element is Xi^/d = —0.5, so that the resonance 
condition J95| ) is satisfied. Arrows on top again indicate zeros of J\. The lattice consists of 20 sites. 

FIG. 20. Magnification of the first avoided quasienergy band crossing seen in Fig. 19. 
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